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The geometry of 4D simplicial quantum gravity with a U(l) gauge field is studied numerically. The phase 
diagram shows a continuous transition when gravity is coupled with a U(l) gauge field. At the critical point mea- 
surements of the curvature distribution of S 4 space shows an inflated geometry with homogeneous and symmetric 
nature. Also, by choosing a 4-simplex and fixing the scalar curvature geometry of the space is measured. 



1. Introduction 

In this talk I am going to discuss on the geomet- 
rical property of space generated by the Monte 
Carlo simulation of 4d quantum gravity mostly 
with one U(l) gauge field. This is the third of 
a series of our reports on this topic at the past 
Lattice conferences |l|,||. Let me first review the 
present status in order to recollect the motivation 
and aim of our work. 

In 1998 there appeared a paper || which influ- 
enced us very much. It is claimed that by adding 
gauge matter fields a new phase, which is called 
the crinkled phase, appears. In some sense this 
result has been already predicted by the analytic 
calculation Q. Taking the conformal mode se- 
riously into account, while the transverse mode 
is treated perturbatively, the string susceptibility 
is shown to have the similar form as the Liou- 
ville theory of 2d quantum gravity. In contrast to 
the 2d case adding matter fields brings the space 
towards a stable configuration. Although the nu- 
merical results gave us a lot of hope for further 
investigation along this line, there still remain a 
few pathological problems. 

At Lattice 98 we added a new data in order 
to clear some of the problems. Firstly, we have 
observed the new phase even in the case with one 



gauge field. Furthermore, we did not have any 
pathological problems previously encountered at 
the transition point. In the Lattice 99 we have 
presented a result of the finite size scaling analysis 
showing that the transition is continuous and the 
fractal dimension measured by two independent 
methods both gave about 4.6 at the transition 
point. 

Although there still remains some points such 
as the phase diagram to be discussed among two 
groups, we feel that in practice we should not 
worry too much about the phase diagrams as far 
as the phase transition remains continuous and 
possibly of second order. 

In the followings we take our numerical results 
seriously, and study the consequence mostly on 
the geometry of space for the case with one U(l) 
field. 

2. Machinery of the numerical simulation 

Let me spend an ample time to describe our 
method employed in the numerical simulation. 
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where c(T) is the symmetric factor of a trian- 
gulation T, and sums are taken over all possible 



triangulations with the topology fixed at S 4 . The 
action is given on the 4d simplicial manifold by 



K4N4 - k 2 N 2 



(2) 



for the Einstein-Hilbert term, where Ni is the 
number of z-th simplex. The two parameters k 2 
and K4 correspond to the inverse of the gravi- 
tational constant and the cosmological constant 
respectively. The action, 
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o(t ijk ){A i3 +A jk + A kl ) 2 



(3) 



corresponds to gauge matter fields, where Aij is 
the non-compact gauge fields defined on the link 
between vertices i and j with the sign conven- 
tion, A^ = —Aji- Sums are taken over all trian- 
gles Ujk in the triangulation T with the weight 
factor o(tjjfe), which is given by the number of 
4-simplices sharing the triangle fcy/.. The term 



AS=-(N 4 -N 4 ) 2 



(4) 



is added to the action artificially in order to con- 
trol N4 around N4 by choosing the parameter 5 
appropriately. 

Measurements are made at every 100 sweeps 
and totally 5 to 10K configurations for each point 
in the graphs. 

3. Results and discussions 

Our primary concern to the present work is 
to make sure that we are dealing with the Ein- 
stein gravity. To this end one may ask oneself the 
following questions: 1) Does the average prop- 
erty of space reflect the classical solution of Ein- 
stein equation? 2) Is there massless graviton act- 
ing between masses? Unfortunately, we do not 
have simulation data with large space configura- 
tion enough to answer the second question, and 
in this report we consider only the first question. 

The first measurement we make is the distri- 
bution of local scalar curvature. The local scalar 
curvature at a 4-simplex s is defined as 



(5) 
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Figure 1. The curvature R S (D) plotted versus 
the geodesic distances D 



where <5 t is the deficit angle of a triangle t, 
S t = 27T-COS- 1 (~ J o(t), 



(6) 



and where Vt and V s are the area and the volume 
of the elementary triangle and the 4-simplex re- 
spectively. In this definition the local curvature is 
an average of curvatures attached to 10 triangles 
t which belong to the 4-simplex s. The lattice 
constant is chosen so that 2Vt/V s = 1. In con- 
trast to the 2d simulation the scalar curvatures 
distribute sharply about average values. 

In order to see how the curvature distributes 
locally in space, we measure the curvature R S (D) 
averaged over a strip between two rings with 
geodesic distances D and D + 1 from a fixed 4- 
simplex s on the S 4 sphere (Fig.l). From the 
figure the space looks homogeneous. Measuring 
R S (D) for several starting simplices we check the 
rotational symmetry of the space. It indeed shows 
that the space is fairly symmetric. 

The Einstain equation with the cosmological 
term has a homogeneous and symmetric solution 
with S 4 topology of the form 



ds 2 = 



dr 2 
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where the parameter tq is related to the scalar 
curvature as R = ■% and to the cosmological con- 

stant as A = if. The space simulated by the dy- 
namical triangulation shows fractal property with 
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Figure 2. The volume V{D) plotted versus the 
geodesic distances D. 

the Hausdorff dimension a little bigger than 4 , 
while the space obtained as a classical solution of 
the Einstain equation is smooth and its dimension 
is apparently 4. 

In order to draw the coordinate frame on such a 
fractal space we need to rescale the distance. By 
comparing the D-step volume, i.e. the number 
of simplices covered within D-steps from a fixed 
simplex, to the volume within a distance r from 
a point of the classical space (eq. ([?])), 

V{r)=^ f-^T- (8) 

J- 



The result is shown in Fig. 2 by rescaling r = 
const x D~ , adjusting the parameters as const = 
0.095 and d H = 4.6. 

The second measurement is for the effect of 
matter. We then check the effect of 'mass' by 
fixing the scalar curvature at a simplex s as 
R s = —1.0. Numerically we reject those geomet- 
rical moves which change R s without any con- 
straint in matter field configurations. The effects 
are shown for the scalar curvature in Fig.l. The 
average property of curvature distribution is sym- 
metric and it can be written as a symmetric so- 
lution of the Einstein equation of S universe: 

(9) 



ds 



1 — r z j(r) 



The scalar curvature of the space is given by 
fl(r) = 12/(r) + 3r/'(r)- (10) 



From this equation it is possible to obtain 
the function f(r) numerically by measuring 
R(r), resulting as 



(11) 



From these numerical results we find that the 
space with and without mass can be considered 
as the solution of Einstein equation with the cos- 
mological term, if we rescale the coordinate ap- 
propriately. However, in order to check whether 
it is the theoretically expected solution of quan- 
tum gravity is left for the analytic calculation for 
the gauge matter energy concentration. 

4. Conclusion 

Taking the numerical simulation of simplicial 
quantum gravity with U(l) gauge fields seriously, 
we have measured the geometry of space. Al- 
though the microscopic structure of space ex- 
hibits fractal, the Hausdorff dimension is 4.6, 
which is rather close to the classical dimension 
4. Furthermore, the scalar curvature distributes 
sharply about its average, and the global geome- 
try seems to have a correspondence to the solution 
of Einstein equation for both cases with and with- 
out mass. These results encourage us to proceed 
further for the correlation measurement in order 
to detect graviton. We also feel that it is now 
worth considering more seriously to make effort 
to compare numerical results to the analytic cal- 
culation along the conformal quantum gravity 
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